Polarized Brillouin Scattering in Salol: Effects of Rotation- Translation Coupling 
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We have studied the 90° VV-polarized Brillouin scattering spectrum of the molecular glassformer 
salol and observed a previously unnoticed VV-dip feature at low frequencies for temperatures above 
300K. This new feature is a consequence of rotation-translation coupling, as recently predicted by 
Pick, Franosch et al [Eur. Phys. J. B 31, 217, 229 (2003)], who showed its relationship to the 
Rytov dip that occurs in the corresponding VH spectrum. The analysis of the spectra shows good 
agreement with the theoretical predictions. 
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Thirty five years ago, the broad quasielastic rotational 
line in the depolarized Iy H (u>) light-scattering spectra of 
liquids composed of strongly anisotropic molecules was 
found to have a previously unknown feature 0,0. With 
high resolution, at temperatures above the melting point, 
a narrow central dip usually designated as the "Rytov 
dip" H is observed. For liquids that can be supercooled, 
with decreasing temperature the Rytov dip disappears, 
the rotational line narrows, and new spectral features de- 
velop symmetrically on the wings and, as the glass tran- 
sition temperature is approached, sharpen to become the 
TA Brillouin components. 

The origin of the Rytov dip was shown by Ander- 
sen and Pecora Q], Q and others || to be rotation- 
translation [RT] coupling, the same phenomenon respon- 
sible for flow birefringence. While the rotational spectral 
line is generated by orientational fluctuations, the dy- 
namics of these fluctuations are modified by coupling of 
orientation to the transverse velocity. 

To explain the evolution of the VH spectrum with 
temperature, the equations of coupled translation-plus- 
rotation dynamics must be modified to include viscoelas- 
ticity, i.e. memory functions must be introduced for 
the relevant transport properties. This synthesis was 
first described in a series of papers by C.H. Wang and 
his coworkers who carried out light-scattering studies of 
a series of molecular glassforming liquids JtJ. If mem- 
ory effects are ignored, Wang's phenomenological theory 
reduces to the Andersen-Pecora results which were ob- 
tained from a two- variable Zwanzig-Mori formalism with 
the Markov approximation. 

Dreyfus, Pick and their coworkers have carried out 
a new phenomenological analysis of the VH spectrum 
which is more complete than previous theories jo]. In 
their most recent publication, this analysis was extended 
to the polarized I V y(u>) spectrum, and it was found that 
this spectrum should also exhibit a characteristic signa- 



ture of RT coupling || . Franosch et al [[[(J implemented 
a parallel Zwanzig-Mori analysis of both the VH and VV 
spectra with results that agree with the phenomenolog- 
ical results of Dreyfus et al and, if memory effects are 
ignored, again recover the VH results of Andersen and 
Pecora. 

In order to test these new theoretical results and, in 
particular, to search for the predicted new effect of RT 
coupling in the VV spectrum, we have carried out an ex- 
tensive light-scattering study of the much-studied molec- 
ular glassforming liquid salol (T M =316K, T G =218K). We 
collected three complete sets of spectra, 90°VV, 90°VH, 
and VH backscattering, at 23 temperatures between 
380K and 210K. A detailed analysis of all the spectra 
will be reported in a forthcoming publication pd]| . In 
this communication we discuss the 90° VV spectra and 
describe the new feature we have observed that confirms 
the predictions of Pick, Franosch et al Q , . In con- 
trast to previous analyses of VV spectra which assumed 
simple addition of the rotational and density-fluctuation 
spectra, RT coupling must be included in the analysis. 

In the phenomenological theory, the conventional equa- 
tions of hydrodynamics are augmented with an addi- 
tional non-hydrodynamic variable to describe the depar- 
ture from isotropy of the average molecular orientation 
within a small volume element. For linear (or axially 
symmetric) molecules, if P(9, <fr, f, t) is the probability of 
finding a molecule at f with its axis u pointing in the di- 
rection (9, (f)), an appropriate variable is the orientational 
density _which forms a symmetric traceless second-rank 
tensor Q: 



,(f,t) = J ddd<f>P(9,<f>,f,t) [u i (9,<l>)u j (9, < f>) - ±8ij] 



(1) 



For the equation of motion of Q, Dreyfus and Pick used 
a damped harmonic oscillator equation 
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QpQij = -urQu - r ® g^Qij + AV ® nj (2) 

where ® denotes convolution, lur is a librational fre- 
quency, r(t) is the orientational friction function, and 
the last term incorporates the coupling of orientation to 
the shear rate Tjj . 

In the Zwanzig-Mori approach, the damped oscillator 
equation for Q is obtained only if both Q and dQ/dt are 
included in the set of selected variables. (With Q alone, 
the equation of motion for Q is purely relaxational.) This 
result parallels the fact that both the density fluctuation 
p q and its time derivative (or current) must be included 
to obtain the usual damped oscillator equation for the 
density-fluctuation (LA) modes. 

The fluctuations of the dielectric tensor include two 
terms, one due to density fluctuations and the other to 
orientational fluctuations: 

Pl{w) — \ up- - bJi - icjr - — 

I Pm 

In Eq. (5), the longitudinal viscosity til(u)) does not 
appear explicitly because we assume that it is propor- 
tional to the shear viscosity rjs(w): t]l(u)) — Vrjs(uj), u>q 
is the bare LA frequency, and Tq is a constant damping 
term due to other processes. 

The RT coupling mechanism enters twice. First in Eq. 
(4) where the first term in the square brackets (S) rep- 
resents the conventional light-scattering channel for LA 
modes via coupling of the dielectric constant to density 
fluctuations (the first term in Eq. (3)), while the sec- 
ond term represents an additional light scattering chan- 
nel provided by coupling of the LA mode to molecular 
orientation. Second in Eq. (5) where the final term gives 
the reduction of the shear viscosity by a factor propor- 
tional to A. For salol, this reduction is ~ 10% and has 
only a minor effect on the shape of the spectrum. 

Experiments: Salol samples for our light-scattering 
experiments were prepared by triple vacuum distilla- 
tion of phenyl salicylate purchased from Sigma Chemical 
Company. The final distillation was made into cylindri- 
cal glass sample cells which were then flame sealed under 
vacuum. The Brillouin scattering experiments were per- 
formed with incident 5145 A monomode laser excitation 
of ~ 250 mW; the spectra were measured with a Sander- 
cock 6-pass tandem Fabry-Perot interferometer. 

In Fig. 1, the upper panel shows the 350K VV spec- 
trum (points) and the VH backscattering spectrum (solid 
line) scaled to overlap the VV spectrum at high frequen- 
cies (30-50GHz). Numerical subtraction of the two spec- 
tra yields the 350K difference spectrum I^tiui) shown by 
the points in the lower panel. For frequencies above ~ 
15GHz, Idif{uj) vanishes as expected. However, for 



6e(f,t)=aS Pm (r,t)I + bQ(r,t) (3) 

Both terms in Eq. (3) as well as their cross products 
must be evaluated to find the complete spectrum. The 
resulting equation for the VV-polarized spectrum Iy V (uS) 
obtained by Pick et al M (with some simplifying approx- 
imations EJ ) is given by 

l vv (u) = ^lm{\R{oj) + —A [S + 1r(lj)] 2 P l (l>)\ 

W [3 p rn 3 J 

(4) 

where the initial f-R(w) term is the pure orientational 
part of the spectrum (which can be determined by 
analysing the VH backscattering spectra), A is the RT 
coupling coefficient, and -Pl(w) is the longitudinal prop- 
agator: 
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FIG. 1. Generating the Idif(^) difference spectra. Up- 
per panel: Iy V spectrum (points) and scaled /yff-back 
spectrum (line) for T = 350 K. Inset: same, with lin- 
ear y and log x scales. Lower panel: the difference 
spectrum obtained by subtracting the scaled /y^-back 
spectrum from the I v ° v spectrum. (The line is I = 0) 



frequencies < 2GHz,I<nf(ui) < 0. As we will see, it is 
this negative region which is the signature of the new RT 
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coupling effect. Idif(ui) is given by Eq. (4) with the pure 
orientation term suppressed: 

I dlf {uj) = -Im(^A [S+%R(u)] 2 P L (w)\ (6) 

UJ I p m 6 ) 

In the fits, we modeled R(uj) = R°Hji((j) and 
Loris{uj) — rf s Hs{uj) with a "hybrid function" 

H x (cj) = [1 - (1 + iwT x )-P* + iujp^t- 1 + iwf- 1 ] (7) 

(where x — R or S) which combines the Cole-Davidson 
function at low frequencies with a tu a power-law at high 
frequencies to approximate the mode-coupling form in 
the frequency region of the experiments |L2| . The param- 
eters in Hr(lo) were determined from the VH backscat- 
tering spectra, and those in Hs{u) and A from the VH- 
90° spectra pJ{ |. Both tr and t$ increase rapidly with 
decreasing temperature with a nearly constant ratio. 

Fits to the density-fluctuation-only model: In 
analyzing the spectra, we initially ignored RT cou- 
pling and fit the difference spectra with a conventional 
density-fluctuation-only model. If RT coupling is ig- 
nored, Eqs. (5) and (6) reduce to 

= -Xm[P L (u)} (8) 

where A = I (q 2 / p m )A.S 2 , and 

P l {lo)=[lo 2 -lu 2 -iluT q -A 2 H s {lo)]~ 1 (9) 

where Hs(uj) is the hybrid function (Eq. [7]) for the shear 
viscosity T]s(l>) = t]° s H s {lu) and A 2 = (q 2 / p m )Vij s . 

We fit the difference spectra to Eqs. (8) and (9), with 
the parameters appearing in H$(io) and A again fixed 
from the fits. This is a conventional analysis of 

Brillouin scattering spectra similar to that carried out 
previously for many other materials (e.g. refs. |]l3|, [H.) 
In the fits, A, A 2 , and uiq were the free parameters. We 
found that good fits could be obtained with fixed Tq 
proportional to T, as was found for propylene carbonate 
jl3| . The choice Tq = T/600 was found to produce values 
of wo quite close to the values computed from ultrasonic 
data. These "density-fluctuation-only" fits are shown in 
Fig. 2 by broken lines. While the fits are generally ex- 
cellent, for temperatures above 280K the low- frequency 
regions of the fits become increasingly less satisfactory 
because the spectra become negative at low frequencies 
while the theoretical fits are always positive. (Note that 
while the difference spectra have low-frequency negative 
regions, 
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FIG. 2. Difference spectra for T =240, 260, 280, 300, 
320, 340, 350, and 370K. The spectra have been scaled 
and shifted vertically for visibility. Solid lines: fits 
to the full theory; broken lines: fits to the density- 
fluctuation-only model. The dotted lines indicate the 
1=0 baseline for each spectrum. Inset: Low-frequency 
region for T = 350K. 

the full Jyv(tj) spectra are everywhere positive). 

Fits to the full theory: Next, we fit the I^n (u>) spec- 
tra using the full theory of Eqs. (5) and (6) with ljq, V, 
and S as the free fitting parameters. The fit results are 
shown in Fig. 2 where we have included both the full fits 
(solid lines) and the density-fluctuation-only fits (broken 
lines). At low temperatures the two fits are indistin- 
guishable, showing that RT coupling is unimportant in 
the frequency range of our spectra. At higher tempera- 
tures, the low-frequency region is fit much better by the 
full theory than by the density-fluctuation-only theory, 
showing that RT coupling, which is ignored in the con- 
ventional density-fluctuation-only analysis, is the source 
of the negative region. 

Referring to Fig. 4 in ref. ||, there is also a tendency 
for the full theory to fall below the density-fluctuation- 
only fits to it at low frequencies when lobt ~ 1. There 
are differences between the lineshapes in that figure and 
ours, however, which we attribute primarily to the use 
of a single relaxation time in their calculations, while for 
salol, tr ~ 15Tg. Because of this large ratio, there is a 
clear separation between the Mountain mode visible as a 
positive feature at temperatures between 260 and 290K 
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and the VV dip, visible for temperatures above 300K. 
The Mountain mode results from coupling of strain to 
structural relaxation and, at low temperatures, is visible 
at frequencies below the Brillouin peak. For T > 300K, 
it is cut off at low frequencies by the VV dip. 

Physically, RT coupling enters because the uniaxial 
strain that characterizes a longitudinal acoustic mode is a 
superposition of simple compression and shear strain. In 
liquids composed of anisotropic molecules, the fluctuat- 
ing longitudinal strain induces a preferential orientation 
of the molecules in the plane perpendicular to q = q'z. 
adding to Se yy which determines the VV scattered in- 
tensity. Referring to the equation of motion for Q (Eq. 
[2]), the steady-state solution gives Qij proportional to 
the strain rate . Therefore, for oscillatory LA waves in 
the low- frequency limit (where lot « 1), if Spij is real, 
then SQij is imaginary. In this limit, the optical coupling 
function C(cu) = [S + (2/3)i?(w)] 2 in Eq. [6] becomes 

[S + ^R(to)} 2 -» [S + iu^RopRTR] 2 « [S 2 + IU^RoPrTr] 

(10) 

and, from Eqs. [5] and [6] in this low-frequency limit, 
I dif (u) - ^[<?V— " S^RcPrtr] (11) 

The second term, being negative, reduces the intensity 
at low frequencies, but it will disappear at frequencies 
lo > Tft where Q can no longer follow the oscillating 
strain. This reduction effect is similar to the Rytov dip 
where, for wtr << l,R (u>) is negative which adds a 
narrow negative component to the center of the broad 
rotational line. 

We therefore conclude that the VV-dip phenomenon 
observed in this experiment is a consequence of RT cou- 
pling as described in the analysis of references H and 
. Finally, we note that there have been many previ- 
ous Brillouin-scattering studies of liquids of anisotropic 
molecules (including salol) , where the VV-dip could have 
been observed. The reason that it has not been reported 
previously is suggested by the upper panel of Fig. 1 
which shows that on a log plot, the low-frequency Ipy(u) 
spectrum is very close to the rotational spectrum. It is 
only when the low-frequency region of the difference spec- 
trum is examined that the small low-frequency VV dip 
becomes apparent. An extensive discussion of this con- 
clusion together with the full underlying data analysis 
can be found in Ref jll]] . 

We thank W. Gotze, T. Franosch, C. Dreyfus, M. 
Fuchs, and C.H. Wang for helpful discussions. This ma- 
terial is based upon work supported by the National Sci- 
ence Foundation under grants DMR-9980370 and DMR- 
0243471. 
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